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The field of quantum noise has been rapidly developing in recent years, 
with the growing possibilities in precision measurements [1], and interest 
in mesoscopic systems as well as in technological applications of physical 
effects at the micrometer and nanometer scales. 

Of particular interest is the study of the statistics of charge transport 
between materials coupled through a contact or through a time dependent 
scatterer. The full statistics of charge transport was studied in a series of 
works by Levitov et al. [2, 3]. This approach yielded interesting results, and 
in particular, they where able to express the full counting statistics in terms 
of a determinant of a single particle operator. Several aspects of Levitov's 
formula (4) where discussed in following papers [4, 5]. 

Our aim in this paper is to present a novel derivation of the original 
Levitov formula. This is done by proving a trace formula (8), which re- 
lates certain traces in Fock space to single particle determinants. Using 
the present approach we find in addition several generalizations, such as a 
corresponding formula for Bosons. 

1. The full counting statistics 

The typical setting is the following: consider particle reservoirs, with given 
temperatures and chemical potentials, which are separated at time zero, 
and are evolving by a second quantized hamiltonian Hq. At some time the 
reservoirs are coupled, through a scattering region, and evolve by a new, 
time dependent hamiltonian. After a time T they are decoupled again and 
one is interested in the statistics of charge transported from side to side, 
i.e. to compute < Q >,< Q 2 > and higher moments. The role of the third 
moment was recently discussed in [6]. 

The full statistics of charge transfer may be conveniently represented by 
the characteristic function of the (charge transport) probability distribution 
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function, denned by: 

X (A l5 ; T) = £ P(<*(t = 0), f3(t = T))e*Xi x ^~^ (1) 

Here the summation is over all states a = (ai, ...), (3 = (/3i, ...) labelling the 
Fock space in the occupation number representation: for fermions these are 
vectors of zeros and ones and for bosons vectors with integer coefficients, 
where «j is the number of particles occupying the single particle state i in 
a. P(a{t = 0),(3(t = T)) is the probability that we started in state a at 
time and finished in a state (3 at time T. Thus terms (a* — /%) appearing 
in the exponent are just the change in the number of particles occupying 
the single particle state i. And the parameters Aj are introduced in the 
standard manner to calculate different moments. By taking derivatives of 
X with respect to Aj, one can calculate arbitrary moments of the charge 
accumulation in state i. For example, 

<Qi>=-ia Ai logxUi,...=o (2) 

And 

(AQ,) 2 =< Qi~ < ft > 2 >= ~dl logxk,...^- (3) 

It is also possible to compute the moments of the charge which is trans- 
ferred to a particular reservoir by taking the derivative with respect to A 
after setting Aj = A for all states % belonging to the desired reservoir. x(A) 
may be interpreted as the reaction of the system to coupling with a classical 
field A which measures the number of electrons on each side. 

For adiabatic change, and short scattering time Levitov et al [2], ob- 
tained the following expression for %: 

X (A) = det(l + n(S^e iqX Se- iqX - I)) (4) 

Where n is the occupation number operator and S is the scattering matrix. 
As was remarked [2], this expression requires careful understanding and 
regularization. In the following we derive this formula in a new manner 
which, we hope will allow a convenient way to address these issues. 
In order to proceed we first write % as a trace in Fock space: 

X(A,T)= (5) 
< a\p \a > | < a\U^\f3 > \ 2 e iq ^i UPi-<*i) 

Where po is the density matrix at the initial time (t = 0) , U is the evolution 
(in Fock space) from time t = to time t = T and a\ , ai are the creation and 
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annihilation operators for a given one particle state i. Here it is assumed 
that the occupation number basis is chosen such that the initial time density 
matrix po is diagonal in it, which implies that the states a are eigenstates 
of the initial Hamiltonian, and measurement of charge in a specific state is 
meaningful. 

Next, we define the second quantized version of a single particle operator 
A (i.e. an operator on the single particle Hilbert space) to be the Fock space 
operator: 

T(A) = ^< i\A\j > a\a r (6) 
Then (5) can be written as: 

X (A,T) = Tr(poe^ Utr(A)U e-^ r W) (7) 

Here A is the matrix diag(Ai, A2, ...). To handle this kind of expressions 
(and to obtain Levitov's formula) we prove in the following section a trace 
formula. 

2. A trace formula 

In this section we prove the following: 

Tr(e r ( A )e r ( B )) = det(l - ^e A e B )^ (8) 

Where £ = 1 for bosons and £ = —1 for fermions (i.e. the creation and 
annihilation operators satisfy Oj-aJ — £aj a j = 

We prove this result for the finite dimensional Hilbert space case, and 
avoid at this point questions regarding the limit of infinite number of states, 
to be addressed elsewhere [7] . 
Proof: 

For an iV dimensional single particle Hilbert space T is a representation of 
the usual Lie algebra of matrices gl(N). Indeed, substituting the definition 
(6), together with the relations obeyed by the creation and annihilation 
operators it is straightforward to check that 

[T(A),T(B)}=n[A,B]) (9) 

is true for bosons and for fermions. By Baker Campbell Hausdorf there 
exists a matrix C such that e A e B = e c . C is an element of gl (N) and is 
given by a series of commutators, since F is a representation, it holds that 



e A e B = e c^ e r { A) e r { B) =e r(C)_ 



(10) 
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Now let us evaluate Tr(e r ^). Any matrix C can be written in a basis in 
which it is of the form diag(/xi, ..fi n ) + K where K is an upper triangular, 
thus we have 

Tr(e r(c 0) = Tr( e r ( dia g(Miv-^n))+r(ft:)) = rp^rxdiagOw, ..,/*„))) = ^-q 
Tr (fj. gwotoi) = _ £ e w)-£ = det(l - ?e c )^ 

(One may also think of T^e 1 "^) as the partition function of a system with 
Hamiltonian — C at temperature fcgT = 1). From this equation (8) follows: 

TV(e r(A) e r ( B) ) = T^rxc)) = det(1 _ ^Cyt = det(1 _ ^AgB)-* (12) 

We remark at this point that this relation can immediately be generalized 
in the same way to products of more then two operators. 
• Let us illuminate our identity with a trivial example: Let 7i be an N - 
dimensional Hilbert space, and choose A = B = 0. Then the dimension of 
the appropriate Fock space is given by 



2 N Femions 
oc Bosons 
as it should be. 



Tr(I) = Tr(e r (°)e r (°)) = det(l - = { 



3. Levitov's formula 

We now turn to give a novel derivation of Levitov's result for the full count- 
ing statistics. In the framework of non interacting fermions the evolution 
U in the expression (5) for \ is j us t the Fock space implementation of the 
single particle evolution U. That means that U f r(A)U = F(U^XU) so that 
by the trace formula (8) for 3 operators, we immediately have 

p -pr(H ) 

x(A , T) = Tr(^-——e^^e-^ r ^) (13) 

Zj 

\ det(l + e- pHo {U^e iqX Ue- iqX )) = 
z 

det(l + n(U^e iqX Ue- iqX - 1)) 

Where Z = det(l + e~^ H °) and n is the occupation number operator 
i+e-^Ho a * the initial time. We note that the result (13) should be viewed as 
the general expression for the counting statistics of noninter acting fermions, 
at any given time, and without any approximation. And may be a good 
start for studying different limits of the problem, as well as regularization 
difficulties. 

Finally, if the scattering time is very small compared with the entire 
evolution, then one may describe the problem in terms of dynamical scat- 
tering operators, S = lim^oo e tHot U(t, —t)e lHot where Hq is the initial free 
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evolution. Using the fact that A commutes with Hq, one obtains in the limit 
of T -» oo: 

X(A) = det(l + n(S^e iqX Se- iqX - 1)) (14) 

Which is Levitov's result (4), as promised. 
We now add a few remarks: 

1. Convergence and regularization: First we note that as long as we 
assume that p exists and has trace 1, then trace of p times a bounded 
operator is also finite, so that x is wen defined. However, problems might 
arise when taking the thermodynamic limit. Taking the infinite volume limit 
may cause the Fock space density matrix to be ill defined (i.e it cannot be 
normalized to trace 1), however, expectation values obtained using it may 
still have meaning. 

If one uses the scattering matrix as in (4) , regularization of the determi- 
nant is needed, since in this case one uses the static scattering matrix as an 
approximation for the true evolution. Indeed, in the limit T — > oo arbitrar- 
ily large charges can pass from side to side 1 , so that the information about 
the length of the time interval has to be put in by hand [2] . The equation 
(13), however, can be shown to be well defined even in the thermodynamic 
limit, for a finite time interval [7]. 

2. Bosons: It is now straightforward to derive an analogous formula for 
the full counting statistics of bosons. The result is simply: 

Xb(A ' T) = det(l-n B (?7te^[/ e -^A_ 1 ))- ( 15 ) 

Where ub is the occupation number operator for bosons. 

3. Rate of charge accumulation: Here we give an example of how one 
may compute the rate of charge accumulation in a box A. We choose a box 
in space, which can be described by a projection P A in the single particle 
Hilbert space (with matrix elements < x\Pa\x' >= S(x — x') if x is in the 
box, and zero otherwise). By setting A = for all states that are outside 
the box A, one finds 

Q A = -id t d x \og X (\,t)\ x=0 = (16) 
qd t Tr(n(U* P A U - P A )) = qTr(n(U^P A U + U ] ' P A U)) = 
qRe < UU ] P A > t 

The angular brackets describe averaging over the distribution at the time of 
measurement t. Charge accumulation is equivalent to current if the box A is 

1 To see this we note that the first moment of (4), which is the transported charge, 
diverges as time goes to infinity if there is a bias between the reservoirs 



6 



I. KLICH 



connected via just one contact to the other reservoirs. This equation should 
be compared to the formula for the current in terms of scattering matrices 
[8, 9, 10, 11] which is of fundamental interest in the field of quantum pumps. 

4. summary 

To conclude, we presented a novel derivation of Levitov's determinant for- 
mula for the full counting statistics of charge transfer. This was done by 
introducing a trace formula (8) which is suitable for translating problems 
of non-interacting particles from Fock space to the single particle Hilbert 
space. The derivation is general enough to allow consideration of new prob- 
lems of counting statistics, in particular, further problems involving bosons, 
or measurement of other operators then charge. We hope that some prop- 
erties of the determinant (13) under various limits, such as adiabatic and 
thermodynamic limits will now be easier to address. 

I am grateful to J. E. Avron for many discussions and remarks, and es- 
pecially to O. Kenneth for his help in the proof of (8). I would also like to 
thank J. Feinberg, L. S. Levitov and M. Reznikov for remarks. 
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